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Abstract. The effective Hamiltonian of strongly correlated electrons on a square 
lattice is replaced by a renormalised Hamiltonian and the factors that renormalise 
the kinetic energy of holes and the Heisenberg spin-spin coupling are calculated 
using a Gutzwiller approximation scheme. The accuracy of this renormalisation 
procedure is tested numerically and found to be qualitatively excellent. Within the 
scheme a resonant valence bond (rvb) wavefunction is found at half-filling to be 
lower in energy than the antiferromagnetic state. If the wavefunction is expressed in 
fermion operators, local SU(2) and U(1) invariance leads to a redundancy in the 
representation. The introduction of holes removes these local invariances and we 
find that a d-wave rvb state is lowest in energy. This state has a superconducting 
order parameter whose amplitude is linear in the density of holes. 


1. Introduction 

Very soon after the discovery of high-T c superconduc- 
tivity Anderson [1] proposed that it was caused by a 
cooperative condensation of carriers moving in a reson- 
ant valence bond (rvb) state of spins. Since then this 
proposal has been studied extensively and the best 
account is in Anderson’s recent lecture notes [2]. Many 
other proposals have been made (for a Review see [3]), 
and there have been questions raised about the use of 
the simplified effective Hamiltonian derived from a 
single-band Hubbard model in the atomic limit that 
forms the starting point of Anderson’s treatment. We 
shall not go into these questions here but just point out 
that two of us (Zhang and Rice [4]) recently gave an 
explicit demonstration that a two-band model describ- 
ing hybridised copper 3d and oxygen 2p states can also 
be reduced to the same effective Hamiltonian in an 
appropriate limit. 

The effective Hamiltonian contains the strict local 
constraint which forbids double occupancy of any site. 
This constraint is very difficult to handle analytically. 
One of the most physically transparent methods to treat 
this type of problem analytically has been the Gutzwil- 
ler approximation which introduces a renormalisation 
of the quantum mechanical expectation values by a 
classical weighting factor [5]. Such renormalisation can 
then be incorporated into a Hamiltonian which may be 
treated by conventional methods. This approach was 
used for the heavy fermion problem by Rice and Ueda 

[6] and was shown to be equivalent to an optimal 
slave-boson formulation by Kotliar and Ruckenstein 

[7] . In this paper we will consider this renormilisation 

t Permanent Address: Institute for Solid State Physics, University of 
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Hamiltonian method for the effective Hamiltonian. 
Unlike the renormalised Anderson Hamiltonian studied 
by Rice and Ueda [6], in the present case the renorma- 
lised Hamiltonian cannot be simply diagonalised and 
we must resort to a further mean-field approximation. 
Mean-field approaches have been considered by many 
authors [8-10], Here we choose to formulate the 
problem in terms of a variational wavefunction. This 
has several advantages. First it shows us that a consis- 
tent mean-field theory must be formulated in terms of 
two expectation values i.e. one must include particle- 
hole amplitudes of the form <c + c> in addition to 
particle-particle amplitudes of the form <c + c + >. This 
point has been recently realised by others as well [1 1]J. 
The coupled equations to minimise the energy have a 
wide class of degenerate solutions at half-filling. Sec- 
ondly a wavefunction formulation is suited to examin- 
ing the role of the redundancy in the fermion 
representation which is not present in the spin represen- 
tation. At half-filling, this redundancy which has its 
origin in the reduction from 4 to 2 degrees of freedom 
per site as one goes from fermion to spin representation, 
appears as a local particle-hole (SU(2)) and gauge (U(l)) 
invariance. The large degeneracy of the mean-field 
description arises from this redundancy and it can be 
shown that it corresponds to the same state in the spin 
representation. Further the appearance of coherence in 
the fermion representation is illusory so that there can 
be no true phase coherence as stressed by Baskaran and 
Anderson [12]. Thirdly, this formulation allows a direct 
comparison with the variational Monte Carlo (vmc) 
results. This allows us on the one hand to check the 

} In [8], Baskaran and co-workers considered the term <c + c> as well. 
However they set it equal to zero in their calculations. 
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validity of the renormalised mean-field theory and on 
the other hand it gives us more insight into the numeri- 
cal vmc results. Both qualitatively and even quantitat- 
ively good agreement is found; for example both point 
to a d-wave paired state as the most stable and a true 
superconductivity order parameter which vanishes at 
half-filling and grows linearly in the deviation from half- 
filling. The largest discrepancy occurs for the anti- 
ferromagnetic state which, within this scheme is higher 
in energy than the rvb state, contrary to the vmc 
results. 

Our treatment is essentially limited to zero tem- 
perature and the extension to finite temperature will be 
non-trivial. Some discussion of the problems of calcu- 
lating excitation energies is given. In particular there are 
two energy scales of excitations given by the gauge 
coherence energy (determined by the kinetic energy) and 
the magnetic coherence energy respectively. Anderson 
[2] has emphasised this splitting of the charged excita- 
tions (holons) and spin excitations (spinons). 

2. The model and the renormalised Hamiltonian 

We study the Hubbard model on a square lattice. In the 
limit of large on-site Coulomb repulsion U and at one- 
half, or, slightly less, filling the Hubbard Hamiltonian 
can be transformed to the form 

H = H t + H s 

H { =-t £ clc M + HC (1) 

<*. J> 1 <* 

h s = j I s t -Sj 

OJ> 

with the local constraint the number of electrons on any 
site < 1. This transformation has a long history, and has 
been used by [13], amongst others. In (1) H t and H s are 
the kinetic and magnetic energies respectively and <i,y> 
represent the nearest-neighbour pairs. S t are the spin = 
\ operators and J = 4 t 2 /U. We neglect terms which are 
higher order in the small parameters t/U and the hole 
concentration S (= 1 — n; where n is the electron 
concentration). 

Since the high-T c superconducting materials show 
strong antiferromagnetic (af) spin correlations [14] it is 
believed that this model contains the essential physics 
for the high-T c superconductivity [2, 4], 

To study the ground state and the excited states of 
(1), we use a projected bcs trial wavefunction as sug- 
gested by Anderson [1] for a rvb state: 

\<P> = P d \<Po> ( 2 ) 

ko>-n (“* + «»* cI t cL fci )| 0 > (3) 

k 

where the Gutzwiller projection operator P d = 
FI, (1 — n ir n n ), and |0> is the vacuum state. u k and v k 
are the variational parameters satisfying the normal- 
isation condition for | <p 0 > : | u k | 2 + | v k \ 2 = 1. 

Some special forms of (2) have recently been studied 
numerically. Using the vmc technique, which treats the 


projection operator exactly [15-19], the energies of 
these states have been numerically calculated. It is 
found that in the square lattice, the projected Fermi 
liquid state (i.e. the state with u k = 0) is unstable 
against d-wave pairing [16]. At half-filling the energy of 
the d-wave state is found [17] to be very close to the 
ground-state energy extrapolated from the exact small 
system calculations [20]. In contrast to the extrapolated 
exact small system calculations, the d-wave trial wave- 
function has no long range antiferromagnetic order [17] 
and may therefore be viewed as an example of a 
quantum spin liquid. A vmc study of superconductivity 
has been made independently by Yokoyama and Shiba 
[19]. They also concluded a possibility of a d-wave 
superconductivity. 

The projected bcs wavefunction is a natural gener- 
alisation of the usual bcs state to strongly correlated 
systems. The projection operator, however, makes diffi- 
culties for an analytic approach. In this paper we shall 
use a renormalised Hamiltonian approach to treat the 
projection operator, and systematically investigate the 
state (2), carrying out explicitly the variational pro- 
cedure. In this approach following Gutzwiller [5] the 
effect of the projection operator on the doubly occupied 
sites is taken into account by a classical statistical 
weighting factor which multiplies the quantum coherent 
result calculated with | cp 0 ). A clear description of the 
method has been given by Vollhardt [21]. The hopping 
energy and the spin-spin correlation of the nearest- 
neighbour sites in the state | <p> are related to those in 
the state | <p 0 ) by 

<cL c j,> = 9,<clc Ja ) 0 

<S, • Sj) = g s (S t • Sj ) o 

where <i4> 0 and < A > are the expectation values in the 
states |<p> 0 and |<p> respectively. The renormalisation 
factors g x and g s are determined by the ratios of the 
probabilities of the corresponding physical processes in 
the states |<p> and |<p 0 >- In figure 1 we illustrate the 
possible hopping processes in these two states. The 
probability of such a process in the state | cp > is 

[”;t(l “ "MtO n jfi U2 

ia) 


j i J i 



Figure 1. The possible hopping processes (a) in the 
non-projected pairing state (3) and (b) in the projected bcs 
state (2). The spins with broken arrows are optional in the 
(a) configurations. 
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while that in the state | (p 0 > is 

CM 1 - «it)«.t( 1 - «jt)] 1/2 - 

n iff are the average electron occupation numbers, (w f = 
I a n io ) which are the same in the states \<p) and | <p 0 >, 
because of the spin symmetry of the wavefunctions. This 
leads to the result [21] 


= 2<5/(l + S). (4a) 

To determine 0 S , we consider the spin exchange 
process shown in figure 2. The spin-spin interaction 
occurs only when both sites are singly occupied. The 
probability for such a process in the state \(p ) is 
(n; T n il n ji n, T ) 1/2 , while in the state | <p 0 > it is 

[n, t (l - n/i)n u (l - - n } t )n it (l - n u )] 1/2 . 

The same result is obtained for the z component inter- 
action S-Sj. Thus one findst 

g s = 4/(1 + <5) 2 . (46) 


It is important to realise that the projection operator 
greatly enhances the spin-spin correlations. To further 
illustrate this point, we list in table 1 all the possible 
two-site states together with their weights and the con- 
tributions to the spin-spin correlation in the half-filled 
case. 

Having determined the renormalisation factors, we 
can define a renormalised Hamiltonian given by 

H' = g { H t + g s H s . (5) 

The energy of the system in the state | <p> can be evalu- 
ated as the expectation value of H' in the state | <p 0 > 

W = (H") 0 . (6) 

Equations (4)-(6) form the basis of our renormalised 
Hamiltonian approach, which is analogous to the 
approach used by Rice and Ueda [6] for the periodic 
Anderson model with the difference that here we make 
a further mean-field approximation. This is because in 
the context of the periodic Anderson Hamiltonian the 
most important physical effect is the renormalisation of 
the f-level to the Fermi surface and not the spin inter- 
action, which would make an exact treatment of the 
effective Hamiltonian impossible. 

To justify this approach, we have carried out Monte 
Carlo calculations. Figures 3-5 show the comparisons 
between the renormalised mean-field theory and the 
essentially exact mc results for these wavefunctions. The 


j 



Figure 2. The spin exchange process in the states (2) and 

(3). 


t In a systematic series expansion on 6 and J/t, the higher-order 
terms of S in (4) should be dropped away to be consistent with the 
effective Hamiltonian (1), where the higher-order terms are not 
included. This, however, does not change the qualitative physics dis- 
cussed in this paper. 


Table 1. This table illustrates the enhancement of the 
spin-spin correlation in the projected bsc state q>, equation 
(2), over that in the bcs state cp 0 , (3), at half-filling. The 
weight of the configurations actually contributing to $ { • S,- 
increases by a factor of four due to the projection. The 
configurations at each site are denoted by 0 (empty state), 
Ti (doubly occupied state), and g (singly occupied state 
with spin g). 


Two-site 





configurations 


Weight 



No of 



Contribution to 

/ i 

configurations 

<Po 

<P 

s;s- or s-s; ? 

0 0 

1 

i 

T6 

0 

— 

G 0 

4 

1 

4 

0 

— 

G G 

2 

1 

5 

1 

2 

— 

G —G 

2 

5 

1 

yes 

ti o 

2 

i 

5 

0 

— 

Ti <7 

4 

1 

3 

0 

— 

Ti Ti 

1 

1 

TS 

0 

— 


quantitative agreement is within 5-15%, while the 
qualitative agreement is excellent for the wavefunction 
( 2 ). 

After replacing the projection operator, the energy 
of the system can be evaluated analytically. The varia- 
tional task is to minimise W in (6). This leads to 
coupled gap equations, which we will derive and solve 
in the following sections. 



Figure 3. A comparison between the renormalised 
mean-field theory (rmf) (see (4)-(6)) and the Monte Carlo 
(mc) result for the kinetic energy <7"> per hole in the 
projected bcs state (2). Both were calculated with a total 
number of sites N s = 82 and a number of holes N h = 8, 16. 
The variational parameter A is related to the parameters of 
the state (2) by 

V k 

u k ~ s k - Ho + [A? + (e k - ti 0 ) 2 T 12 

where g 0 is a parameter, and e k is given by (8b). In the 
d-wave pairing state, A k = A (cos k x — cos k y ) and in the 
s-wave state, A k = A. The full circles and squares are the 
mc results for the s- and d-waves respectively. The dotted 
and broken curves through the mc results are guides for the 
eyes. The second pair of dotted and broken curves are the 
results from rmf for s- and d-waves respectively. 
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Figure 4. A comparison between the renormaiised 
mean-field theory ((4)-{6)) (rmf) and Monte Carlo (mc) 
results for the nearest neighbour spin-spin correlation 
<S, ■ S y > in the projected bcs state (2). Both were calculated 
with a total number of sites A/ s = 82 and a number of holes 
A/ h = 0, 8, 16. The variational parameter A is related to the 
parameters at the state (2) by 

v* A* 

u k ” s* - Ho + [A* + (e k - n 0 ) 2 y 12 ' 

where pi Q is a parameter and e k is given by (86). In the 
d-wave pairing state, A k = A(cos k x - cos k y ) and in the 
s-wave state A k = A. The full circles and squares are the mc 
results for the s- and d-waves respectively. The dotted and 
broken curves through them are guides for the eyes. The 
second pair of dotted and broken curves are the results 
from rmf for the s- and d-waves respectively. 


3. Gap equations 

In this section we derive the gap equations for the pro- 
jected bcs wavefunction within the renormalised Hamil- 
tonian scheme described in §2. We consider only the 
even-parity case, i.e., u_ k v± k = u k v f, and |u_ k | 2 = 
I 2 . 

Evaluating (6), we obtain 

W = 2 g t Ze k \v k \ 2 + N s -> 

k 


x Z v k-k( 1 v k I 2 1 v k , I 2 + u k vtv k . «f.) 

(7) 

k,k' 


where N s is the total number of sites, and 


•« 

cn 

1 

II 

(8fl) 

«*= ~t7k 

(8b) 

y k = 2(cos(/c x ) + cos(k y )) 

(8c) 



Figure 5. The nearest neighbour spin-spin correlation 
function <5, • S y > as a function of electron filling in the 
projected Fermi liquid state (u k v k - 0 in (2)). The 
renormaiised mean-field theory ((4)-(6)) (rmf), broken 
curve, agrees well with the Monte Carlo (mc), full circles, 
result in the entire filling region. 


Note that e k and V k have the same functional form, 
since ff s is derived by kinetic exchange. The electron 
number operator N = L k(T c k(r c ka has expectation value 
<iV> = 2E k | t; k | 2 . 

Let ^ be the chemical potential of the system, the 
quantity we want to minimise is 

W= <H' - fiN}. o . 

While minimising W with respect to u k and v k for fixed 
ju, one must remember that H f is also a function of 3 
due to the renormalisation factors, and hence a function 
also of v k . Carrying out this procedure, we find that 



l«*l 2 = i(l + £*/£*) 

(9a) 


*• 

II 

i 

*Tn 

XT 

bn 

(9b) 

and 

u k v k ~ 

(9c) 

where 

£* = «* 2 + l^| 2 ) 1/2 . 

(9 d) 

The parameters A k and £ k are dimensionless, and 
they are related to the particle-particle and particle- 
hole pairing amplitudes respectively 


& k = A x cos (k x ) + X, cos (k y ) 

& = h-Z x cos (k x ) - i y cos (k y ) 

(10) 

where 

h = (g t e k ~ Mh.J) 


and p is related to the chemical potential \x by 
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In (10) 

= ( C ?T C J+ti — c UtCj+ T f) 0 (12) 

t' = Z<ct.c l+ ',.> 0 (13) 

<7 

with t = x and y, i + r denotes the nn of i in the t 

direction. Since we consider the even-parity case, is 

real, but X T can be complex. A k and satisfy the follow- 
ing coupled gap equations ; 

A* - AT.- 1 1 y*-*.a*'/(2 £*•) 

fc' 

Z k = h + N s - l Zy k - k 'tv/(2 E k ). (14) 

k' 

The first one is the same as the usual bcs gap equation. 
The second one originates from the particle-hole correl- 
ation. From (12), it is clear that A k is related to the 
pairing in the unprojected state | <p 0 }. It describes the 
‘smearing’ of the pseudo-Fermi surface. However, A k is 
not the superconducting order parameter in the project- 
ed state | <p> in our theory. E k turns out to be the 
quasi-particle excitation energy (in units of %g s J) in the 
pairing state as we will show in §5. 

The coupled gap equations (14) are the basic equa- 
tions in our approach. They can also be written in the x 
and y component form : 

A z = /V s ~ 1 Z cos(k t )&JE k 

k 

-N; 1 I CO* WE*. (15) 

k 

The gap equations must be solved simultaneously with 
the hole concentration equation, S = £* £ k /E k . 

Before we discuss the non-trivial solutions, we note 
that A k = 0 is a trivial solution of the gap equations. 
This corresponds to the projected Fermi-liquid state. In 
this case, £ k changes sign at the Fermi surface. The 
parameters £ x = £ y ( = £) are given by 

f-Ar." 1 I cos (k x ) + cos (k y ). 

The volume of the Fermi sea is determined by the 
number of electrons. The energy per site is 

w = -4 g,t£ - %g s JS 2 . 

In particular, w = — 48 /n*J zz — 0.49 J in the half-filled 
case. It will be shown in the next section that this trivial 
solution is unstable against the pairing states with 

A k # 0. 


4. Solutions of the gap equations — half-filled case 

At the half-filling, <5 = 0, ]x = 0, and there is no kinetic 
energy. We are interested in the possible lowest energy 
states of the solution of (14). The total energy of the 
system has a simple form in this case by use of (7) and 
(14) 

w = -f g s JZE k . 


Therefore the lowest energy states correspond to the 
maximum value of E k E k . For this reason, we use an 
ansatz for E k to examine the solutions of the gap 
equations! 

E k = C(cos 2 (k JC ) + cos 2 (ky )) 11 2 (16) 

where C is a parameter to be determined. Note that 
such a choice gives only four point zeros for E k . The 
gap equations then reduce to a single equation, and we 
get 

c = Z (cos 2 (fe,) 4- cos 2 (k y )) 112 

which has numerical value C ^ 0.479. The energy per 
site is 

w = -$g s JC 2 = - 0.688 J. 

This energy is much lower (about 20%) than that of the 
projected Fermi-liquid state found above. The para- 
meter describes the pairing correlation in the renor- 
malised Hamiltonian. Finite values of A k indicate the 
binding of the electron pairs in the pairing states (3). 

We now determine the parameters \ and £ T 
required for the choice of E k in (16). Using (9 d) and (10), 
we find that they should satisfy the following simulta- 
neous equations: 

i 2 x + \K \ 2 = $ 2 y+\K \ 2 = C 2 

2^ y + (M? + HC) = 0. (17) 

There is a wide class of parameters which satisfy the 
conditions (17). All the states in this class give the same 
expectation value of the renormalised Hamiltonian H' 
in (5). Therefore at the half-filling, H f has a large degen- 
eracy of ground states. For real A x and A y , these states 
can be illustrated diagrammatically as shown in figure 
6(a). A few examples of these states are 
d-wave pairing : 



= i, = C/y/ 2 

(18a) 

d-wave density matrix : 

a,-*,-*.- 

-Zy = C/j2 

(18*) 

chiral state : 

A x = —i& y = C 

o 

II 

p*. 

II 

X 

(18c) 

anisotropic state: 

K = $y = C 

a, = fx = o- 

(18<0 


We remark that the d-wave density matrix state is dif- 
ferent from the extended s-wave state proposed by Bas- 
karan, Zou and Anderson [8]. In their theory, the 
particle-hole amplitude <c + c> 0 is not included, i.e. 
£ x = = 0. Therefore their state has the same energy 

as the projected Fermi-liquid state. The d-wave pairing 
state was studied numerically in [17] and [18], and the 
chiral state was discussed in [10]. They belong to the 
solutions of the same gap equations in the present 
approach. 

+ There might be other solutions for the gap equations. We believe 
that the form of the energy given by ( 16 ) gives the lowest energy. 
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(a) (ft) 

$ 



Figure 6. (a) Diagrammatic illustration of the degenerate 
ground states for the renormalised Hamiltonian (5) at 
half-filling. ^ A = & x x + & Y p, with f t , S t given 

by (12) and (13). The full arrows represent the d-wave 
pairing state with % 1 A, and |£| = | A| . All the states in (17) 
with the real parameters A t may be obtained by rotating £ 
and A simultaneously by angles 6 , or reflecting the two 
vectors about the i axis, (b) An illustration of the SU(2) 
gauge transformation. A state described by (£', A') in (a) 
can be obtained by a local SU(2) from the state (£, A), under 
which at the four sites of the plaquette transform 
according to (19), with (a, , /?,) as denoted. The minus sign in 
front of the parentheses corresponds to the states in (a) 
after a reflection about the* axis. The transformation 
operator c io at other lattice sites is determined by a 
translation. 

The degeneracy of the ground states of H’ may be 
explained using the local SU(2) symmetry of the Heisen- 
berg Hamiltonian H s . This symmetry has been pointed 
out by Anderson [2]. Very recently Affleck and 
co-workers [22] have studied the invariance under a 
time-dependent SU(2) gauge transformation, and dis- 
cussed the equivalence of some different mean-field 
theories. Here we wish to show that the symmetry 
described in (17) is a sub-group of the local SU(2), 
which preserves translational invariance and even parity 
that we imposed in deriving the gap equations. 

Consider a local SU(2) gauge transformation, under 
which electron operators at site i transform as 

c h a i C ?r + Pi c ii 

c n -* ~P* c h + a * c n ( 19 ) 

where a t af + = 1. These are the particle-hole 

transformations with spin conservation. H s , hence H f at 
half-filling is invariant under these transformations. 
Therefore all the states related by (19) are degenerate. A 
set of the parameters X t and in our gap equations 
transforms to another set of the parameters under (19). 
The transformations corresponding to the degenerate 
state (17) derived from the gap equations are represent- 
ed in figure 6(b ) , for real values of £ r . There is also one 
to one correspondence to transformations of the d-wave 
pairing state to a state with complex A t such as the 
chiral state (18c) with a and /? complex. 

Local U(l) gauge symmetry is a sub-group of the 
local SU(2) symmetry since it is of the form 

Ci* c i9 = c f<T exp(i Ot). 

A general choice of 6 t transforms the bcs pairing state 
(3) to a non-Bes-type state, which has the same energy 


as the former. Since under such a transformation the 
Bloch coherence is lost, i.e. a state of the form 

C k , ' = £ exp(i£ • Rj)c h „ = £ exp(i* • Rj + i 0)c u „ 
j J 

is no longer a coherent superposition of the original 
states. Yet we can equally well pair c* t with cl fci , and 
the energy would be the same as if we pair c* T with 
ct kl . Thus the apparent coherent fc-space pairing in the 
bcs wavefunction (3) is illusory. The absence of a coher- 
ent pairing order parameter at half-filling as a conse- 
quence of the U(l) gauge invariance has been stressed 
by Baskaran and Anderson [12]. 

It is important to realise that the states that are 
degenerate due to the SU(2) gauge symmetry are the 
unprojected states | <p 0 > of (3), rather than the physical 
states | <p > which obey the strict local constraint. How 
does a projected state change under SU(2)? According 
to (19), a vacuum state (empty state) at site i must trans- 
form under SU(2) as 

1 0>, -» — ft? cf T cff 1 0>( . 

This ensures that the vanishing of the state c,J0>, 
remains unchanged under the transformation as 
required physically. However a singly occupied electron 
state transforms under SU(2) as 

1 0>( -* e w ‘cj„ I O), . 

At half-filling, each site is singly occupied. Therefore 
any half-filled state | <p> transforms into itself except for 
an overall phase factor under the SU(2) operator U : 

i 

Although U does not commute with the projection 
operator P d , we observe for the half-filled state | <p 0 > 

up d ]<p 0 > = p d Oi<p 0 >. 

One way to see this is to notice that there is no empty 
site state also in P d 1 4> 0 ). Thus we can rewrite 

p a \<Po > = n ("it - »ii) 2 i<Po>- 

i 

The SU(2) transformations all commute with the 
operator (n iT - n H ) 2 . Let | <p} = P a \<f> 0 }, and |^' 0 > = 
U | <p 0 ), then 

Pd | cp'o> = p d 0 1 <^o> = 0 1 <P> = e* | = e ie P d | cp 0 y. 

This proves explicitly that the two states | <p 0 > and | <p' 0 ) 
related by any local SU(2) gauge transformation corre- 
spond to the same state | <p) except for a phase factor. 
Therefore all the states in (17) correspond to the same 
physical state. The RVB ground state is non-degradable. 

There are also redundancies in the higher energy 
states in the fermion representation. For instance, the 
state of Baskaran, Zou and Anderson [8] at half-filled is 
identical to the projected Fermi liquid state, because the 
former transforms to the latter under (19) with a = 1, 
P = 0 in one sublattice, and a = 0, /? = 1 in the other. 
This equivalence was also pointed out by Yokoyama 
and Shiba [18] in a different way. 
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We now comment on the local gauge symmetry in 
the original Hubbard model, which in terms of the orig- 
inal fermion operator is 

H H = -t X (4 d ja + HC) + u X 4 t 4 t 4, 4 • 

<» 7 >. * 

This Hamiltonian is not invariant under local gauge 
transformations with respect to the operators d ia . 
However, up to any finite order in t/U there exists a 
canonical transformation, which eliminates the doubly 
occupied sites [13] 

// H -H eff = e iS H H e- i5 . 

At half-filling, all the odd order terms in t/U vanish. 
The Hamiltonian is locally gauge invariant with respect 
to the electron operators in the new representation, i.e. 
the c ia of (1) are Wannier operators of the old represen- 
tation, d ig . Therefore the local gauge symmetry holds to 
any finite order in perturbation theory in t/U . This is 
the same as saying that the system has undergone a 
transition to a Mott insulator [13]. 

We have so far only examined the projected 
Bcs-type trial wavefunctions. It is likely that the true 
ground state of the model Hamiltonian (1) at half-filled 
is the af state [19]. Recently, Yokoyama and Shiba [18, 
19] have studied a projected Hartree-Fock-type af 
state. Using vmc they found the energy per site at 
half-filled to be — 0.642 J, slightly lower than — 0.636 J, 
the value found by Gros [17] in the d-wave pairing 
state by using a similar technique. However we may 
argue that the holes favour the pairing state away from 
half-filled because of the gain in kinetic energy. We have 
also applied the renormalised Hamiltonian approach to 
the af states. Within this approximation, we find that at 
the half-filled, the af state has higher energy than the 
rvb state (2), in contrast with the vmc results. We 
present the derivation and the results in Appendix 1. 


5. Non-half-fllled case 
5.1. Ground state 

Firstly we examine the energy needed to create propa- 
gating Bloch states. The simplest states for holes have 
the form 

\^> = c ia \cp) 

which destroys a real electron at site i . We may also 
make a propagating Bloch state for the hole of the form 

\4> P <r>=c pg \<py (20) 

A rigorous calculation is possible at half-filling. 

We consider any translationally invariant spin 
singlet state | <p> at half-filling. Let us denote by a the 
nn spin-spin correlation in | cp >, a = <S, * Sj > <p . Then 
the magnetic energy loss of a hole in the state <£ ifr is 
— 4a J, because the four bonds connecting the site i are 
mixing. Since the matrix of H s in (2) between any states 
|d> i<y > and | <!>;*> is diagonal, the moving hole state of 
(20) has the same magnetic energy as in | 


The kinetic energy of the hole in (20) is given by 

<H,) P = 2 f £ <n ia n j<7 + S^S~) V exp(i p ■ R j{ ) + hc 

<y> 

where < denotes the expectation value in the half- 
filled state. Using the fact that = 1, we get 

<H t > p = t( 1 + 4a)(cos (p x ) + cos (p y )). (21) 

Since a -0.33 for the ground state, (21) gives a band 
width for a Bloch hole of 0.64 1 1 \ . The minimum energy 
to remove an electron and create such a Bloch hole is 
— 0.32 t — 4a J. 

The Bloch states are not the lowest energy states of 
the holes however. We now apply the gap equations to 
study a system with a few pair of holes. The energy to 
create a pair of holes is —2fi by the definition of 
the chemical potential. Since the parameter Ji = 0 at 
the half-filled, (11) gives the energy per hole to be 
N~ l (dH'/dd) 0 , a quantity related to the unprojected 
state |<Po> at half-filling. In the presence of holes, the 
kinetic part of the Hamiltonian explicitly breaks the 
SU(2) and U(l) gauge symmetries, while the Heisenberg 
spin part remains invariant under these symmetries. 

Using (4) and (5), for the states | <p 0 > described by 
(17), the magnetic energy per hole is 



a value equivalent to the loss of four bonds in the spin- 
spin correlations, and it is the same for all these states 
as a consequence of the SU(2) gauge invariance of the 
spin part of the Hamiltonian H s . The kinetic energy per 
hole in this case is given by 

T= w. (i)<"'>° 

is the particle-hole correlation in | <p 0 > as defined in 
(13). When the holes are introduced, a fraction (g/) of 
this correlation becomes coherent in the state | <p>. 
Therefore the larger values of correspond to the 
lower kinetic energy of the holes. But are subject to 
(17). The kinetic energy can be written in the following 
form by using (17): 

T = — 4t(2 C 2 - \A X + | 2 ) 1/2 sgn(^ + y. 

Different parameters A t and c t describe different states 
with different energies upon doping. The above hole 
kinetic energy expression immediately leads to the 
important conclusion that the d-wave pairing state, 
where = 0 gives the best kinetic energy, which 

is 

7iwave= -4(v/2)tt= -2.7U. 

Both the d-wave density matrix state (18 b) and the 
chiral state (18c) have zero kinetic energies, and are not 
favoured upon doping. The kinetic energy for the 
d-wave pairing state in our analytical approach is quite 
close to the vmc result, where it is found to be —2.55 t 
for systems with 10% holes [17]. It is also substantially 
below the value found for a Bloch hole. 
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The introduction of some holes breaks the local 
gauge symmetry and causes the ground state to be 
coherent. The stable lowest energy state upon doping is 
the d-wave pairing state. 

It is worthwhile remarking that in a Hubbard model 
it is the expectation value of the kinetic energy which 
determines the integrated optical weight associated with 
the charge carriers in the /-sum rule [23] and in the 
present case this optical weight is proportional to the 
number of holes with an optical mass determined from 
the proportionality constant of order at ~ 1 (where a is 
the lattice constant). 

The gap equations for the finite hole concentrations 
can be solved numerically. Here we shall consider only 
the most stable state — the d-wave pairing state. In this 
case, we set A x = — A y = A, and £ x = = f. The four 

equations (13) reduce to two because of the symmetry 
between the x and y components. These two equations 
uniquely determine A and £ for the fixed values of t , J 
and /2. The numerical results of the gap equations are 
plotted in figure 7 for A as a function of the hole con- 
centration. 

We now discuss the superconducting order param- 
eter. As mentioned in §3, A is not the order parameter. 
The superconducting order parameter is 

As c(Rij) = ( c h c ji ~ C U c )\ ) 

an expectation value in the projected state (2). This 
quantity describes the Cooper pairing in a real space 
representation. We shall adopt the Gutzwiller method 
to calculate this quantity. In analogy to the derivation 
for the hopping energy in §2 we find that the nearest- 
neighbour sites i and j 

( c h c Ji) = 9ti c h c Ji>o • 

Therefore for nearest-neighbour sites, the order para- 
meter is related to the variational parameter A in the 



Figure 7. Variational parameter A and superconducting 
order parameter A sc as functions of the hole concentration 
S for a choice of t/J — 5 in the d-wave pairing state. 


gap equations by 

As c = 9t&- (22) 

The value of A sc as a function of 5 is plotted in figure 7 
in comparison with A. A sc vanishes linearly near <5 = 0. 
A sc found in our theory is in good agreement with the 
Monte Carlo results [17]. The absence of the supercon- 
ducting order parameter at half-filled obtained from (22) 
agrees with the discussion in §4 from the viewpoint of 
the local gauge symmetry. 

The kinetic energy of holes in the a f state is found 
to be quite high in our analytic approach (see Appendix 
1). However, the Gutzwiller approximation we adopted 
is too rough to determine whether a rvb or a f state has 
lower energy. Numerical results of vmc [17-19] suggest 
both the spin-spin correlation energy and the kinetic 
energy of the holes between the d-wave pairing state 
and the af states are very close. The question which 
state is more favourable in energy remains unresolved. 

5.2. Excited states and finite temperatures 

We begin by examining the spin degrees of freedom in 
half-filled and near half-filled cases. An excited state can 
be created by applying the spin raising operator to a 
specific site to obtain 

|4' j , + > = S?P d |<p 0 >. 

We can commute Sj with P d to obtain 

l'¥ i ' + > = P i S!\<p 0 > 

= £ exp[i(p - p’) • P,)]P d c ; t cyj <j9 0 >. 

P- P' 

This state is therefore a superposition of two indepen- 
dent quasi-particle states similar to a metal where the 
low energy excitations are made up of superpositions of 
electron and hole states. The quasi-particle states 
(spinons) can be defined by 

1^) = ^ FI (u k + v k el t cl kl m. (23) 

k*p 

The quasi-particie energy E p is defined to be the differ- 
ence of the expectation values of K = H - pN in the 
state liA pt > and in the ground state j <p). We use the 
Gutzwiller method to calculate the energy of the state 
(23). The energy difference between the two states con- 
tains two parts. One is due to the changes of the renor- 
malization factors g t and g s , the other comes from the 
change of the wavefunction itself. The former just 
cancels exactly the second term in p in (11). Using (7) to 
calculate the energy difference due to the wavefunction 
change, we get 

E P = (1 - 2t s p + N^ £ V k . p v 2 k - 

+ 2u k v k N s -^V k . p utv k . 

k 

Applying the gap equations to simplify the expression, 
we obtain 

E p = $g s J E p . (24) 
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Note this energy is independent of the local SU(2) gauge 
and does not depend on the particular fermion repre- 
sentation. At the pseudo-Fermi surface, where by defini- 
tion = 0, we have E p = \A p \. Since the state (23) 
breaks a pair of electrons, E p describes the binding 
energy of the pair at the pseudo-Fermi surface. The 
excitation energy depends on the particular rvb. For 
the Fermi liquid state ( u k v k = 0) then E p = 0 over the 
whole pseudo-Fermi surface. However in the ground 
rvb state it vanishes only at four points, e.g. when 
n = 1, E p = E p = 0, if (p x , p y ) = ( ± jt/2, ±n/2) and the 
density of spinon states at low energies is 

2E 

We turn now to a brief discussion of the system at 
finite temperature. The extension of the mean-field gap 
equations to finite T is not so straightforward. The exis- 
tence of a finite A k is controlled by the energy scale of 
E k i.e. by J in (24). On the other hand if we consider the 
limit <5 ^ 1 there is a very small energy scale oc<5r which 
controls the definition of a coherent gauge. In other 
words it is only the kinetic energy which allows us to 
determine the gauge uniquely and at temperatures 
J > T > 8t, the gauge coherence will be lost. Yet in this 
temperature range the magnetic coherence survives 
since as we have stressed earlier this is independent of 
the choice of gauge on each site. The properties of the 
system in this temperature region are clearly very differ- 
ent from Fermi liquid behaviour as Anderson has 
stressed and these two energy scales should correspond 
to his ‘holon’ and ‘spinon’ energy scales respectively. 
The thermopower should obey the Heikes formula [24] 
and we can expect only a low mobility of the holes. 
However, a more detailed study of this regime is 
required. 


6. Discussion 

We have used a variational method to study a projected 
bcs trial wavefunction for the square lattice effective 
Hamiltonian. Using a Gutzwiller approximation to 
treat the effect of the projection operator, we obtained a 
renormalised Hamiltonian in which the projection oper- 
ator is replaced by renormalisation factors. This 
approximation is shown to be in good agreement with 
numerical Monte Carlo calculations for such projected 
wavefunctions. In this mean-field approximation both 
particle-particle and particle-hole pairing amplitudes 
must be included. The fermion representation for the 
ground state at the half-filled band is highly redundant, 
due to a local SU(2) invariance at exactly half-filling. 
This redundancy is reflected in an apparent degeneracy 
of the bcs trial wavefunction before projection. Doping 
destroys the local SU(2) invariance and splits these 
degenerate states, and we find that the stable state upon 
doping is the d-wave pairing rvb state. In this rvb state, 
electrons are paired even at half-filling and it costs an 


energy of order J to break a pair. These pre-existing 
electron pairs lead to a non-zero superconductivity 
amplitude upon doping, and the magnitude of this 
superconducting amplitude or order parameter is 
shown to be proportional to the hole concentration 8 
when 8 is small. The elementary excitations at 
half-filling are the projected bcs quasi-particle states or 
spinons, with four point zeros on the pseudo-Fermi 
surface. 

Our analytic approach can also be applied to id 
and systems with dimensionality d ^ 3. We find that 
lowest energy state in id is the projected Fermi liquid 
rvb without electron pairing, as shown in Appendix 2. 
Our theory predicts no superconductivity in a ID rvb. 
For large d , the energy per bond in the rvb pairing state 
is proportional to 1/d , reduces relative to an af. So the 
pairing state is particularly favourable in 2D. The precise 
form of the 2d phase diagrams which depends sensiti- 
vely on the relative energies of the af and d-wave rvb 
states as a function of 8 is too subtle a question to be 
settled by the approximation we use here. 

There are many questions that require further inves- 
tigation such as the exact relationship between the dis- 
cussion here in terms of phase coherence among the 
paired electrons and Anderson’s ‘holon’ [25] concept or 
the nature of the high-temperature phase where this 
phase coherence is lost but strong magnetic correlations 
remain and presumably do not lead to a Fermi liquid 
that is the usual description of a normal state. 
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Appendix 1 : The projected spin-density-wave state 

In this Appendix, we use the renormalised Hamiltonian 
approach to study the projected spin-density-wave state 
for effective Hamiltonian (1). That state was proposed 
and studied using vmc by Yokoyama and Shiba [18]. 
The generalisation of the Gutzwiller method to the anti- 
ferromagnetic states for the hopping process was formu- 
lated by Ogawa and co-workers [26]. 

The projected spin-density-wave state [18, 19] is 

l«A> = ^dl^o) (Al.l) 

I *Ao> = n ( u k c L + ™*cJ + g,<t)|0> (A 1.2) 

ki 7 

where k runs over the Fermi sea, Q = n/a (1, 1), and 
«* = [(! + cos 9 k )/2y 2 
v k = [(1 - cos 0*)/2] 1/2 
cos 9 k = y k /( A af + y*) 1/2 . 

A af is a variational parameter, and y k is given by (8c). 
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In a study of the expectation value in the state 
(A 1.1), we use the Gutzwiller approximation to replace 
the projection operator by renormalisation factors. In 
analogy to the analysis we discussed in §2, we find that 

1 -n 

1 — 2n t njn 
9, = (1 — 2n t njri)~ 2 

where n t and n l are the spin-up and spin-down electron 
occupation number of state |^ 0 ) i n one sublattice 
respectively. The renormalisation factors reduce to (4) in 
the case n T = n A , and the form for g t agrees with [26]. 

Within this scheme, we obtain the energy per site at 
half-filling 

w = — 2J(a 2 + 6b 2 )/(l + a 2 ) 2 

where 

a = N 7 1 Z a af/(Aaf + yl) 112 

k 

b = (SN a r i ly 2 M F + yl) 112 - 

k 

The spin-spin correlation <5 ; • S)> = ^w, and the stag- 
gered magnetisation is 

Af s = n T — n t = 2a/(l + a 2 ) 

with h a the occupation number in the state It/'). 
<5j • Sj) and M s are plotted in figure A1 as functions of 
A af . Figure A1 also shows the vmc calculations [18], 
The case A af = 0 corresponds to the projected Fermi 
liquid state, while A af -* co corresponds to the Neel 
state. The results agree well for large values of A af , but 
there are substantial deviations for small A af . 

The kinetic energy per hole in our theory is 

T = — 16tb/(l +a 2 ). 

For the optimal value of A af ( ~ 0.9), T — -2.16 1, sub- 
stantially higher than that in the d-wave pairing state. 
Note that this value is also higher than that found in 



Figure A1. Spin-spin correlation and staggered 
magnetisation M $ as functions of A af in the projected 
spin-density-wave state. The full curves are the results of 
the renormalised Hamiltonian approach, and the broken 
curves are the vmc results (extrapolated to the infinite 
systems) by Yokoyama and Shiba [18, 19]. 


the vmc calculation [18], where the optimal A af is 
found to be much smaller. 


Appendix 2. rvb in a id system 

The renormalised Hamiltonian approach can be 
straightforwardly applied to the model (1) in id. Using 
the projected bcs wavefunction (2), and the same tech- 
nique for 2D, we have found that the rvb ground state 
at the half-filling is described by an equation between £ x 
and A x (defined in ( 1 2)— ( 1 3)) : 

I K\ 2 + £ = C\ (A2.1) 

with 

C 1 =(21V s )- 1 Xlcos(k x )| = 2/7t 2 

kx 

(A2.1) is parallel to (17) in 2D. Similar to the 2D case, 
different parameters in (A2.1) are related to each other 
under the SU(2) gauge transformation, and correspond 
to the same physical state. This state is described by the 
projected Fermi liquid state, where A x = 0, = C x . 

Earlier Monte Carlo calculations [27] and more recent 
exact calculation [28] with this wavefunction have 
shown that the energy of this state is extremely close to 
the exact solution [29]. The energy per site in our 
analytic mean-field approach is — 6/n 2 . This value devi- 
ates by about 37% from the true result [27, 28]. This 
quantitative discrepancy however is not surprising, 
because the Gutzwiller approximation is poor in id. 

In parallel to the discussion in §4, we can study the 
system with some holes. We found that the stable 
lowest-energy state corresponds to = 0. Introducing 
the finite value of A x , the system loses kinetic energy. 
Therefore we expect there is no electron pairing and no 
superconductivity in this id rvb. 
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